SCIENTIFIC 

REPORTS 




SUBJECT AREAS: 

SPINTRONICS 

MAGNETIC PROPERTIES AND 
MATERIALS 

FERROMAGNETISM 

ELECTRONIC PROPERTIES AND 
MATERIALS 



High domain wall velocities via spin 
transfer torque using vertical current 
injection 

Peter J. Metaxas 1,2 , Joao Sampaio', Andre Chanthbouala 1 , Rie Matsumoto 1,3 , Abdelmadjid Anane 1 , 
Albert Fert 1 , Konstantin A. Zvezdin 456 , Kay Yakushi ji 3 , Hitoshi Kubota 3 , Akio Fukushima 3 , Shinji Yuasa 3 , 
Kazumasa Nishimura 7 , Yoshinori Nagamine 7 , Hiroki Maehara 7 , Koji Tsunekawa 7 , Vincent Cros 1 
& Julie Grollier' 



Received 
8 November 2012 

Accepted 
25 April 2013 

Published 
14 May 2013 



Correspondence and 
requests for materials 
should be addressed to 
P.J.M. (peter. 
metaxas@uwa.edu. 



1 Unite Mixte de Physique CNRS/Thales and Universite Paris-Sud 11,1 Ave. A. Fresnel, 91 767 Palaiseau, France, 2 School of 
Physics, M01 3, University of Western Australia, 35 Stirling Highway, Crawley WA 6009, Australia, 3 National Institute of 
Advanced Industrial Science and Technology (AIST) 1-1 -1 Umezono, Tsukuba, Ibaraki 305-8568, Japan, 4 A. M. Prokhorov General 
Physics Institute of RAS, Vavilova str. 38, 1 1 9991 Moscow, Russia, 5 Moscow Institute of Physics and Technology, 141 700, 
Institutskii, 9, Dolgoprudny, Moscow Region, Russia, 6 lstituto P. M. s.r.l., via Cernaia 24, 1 01 22 Torino, Italy, 7 Process Development 
Center, Canon ANELVA Corporation, Kurigi 2-5-1 , Asao, Kawasaki, Kanagawa 2 1 5-8550, Japan. 

Domain walls, nanoscale transition regions separating oppositely oriented ferromagnetic domains, have 
significant promise for use in spintronic devices for data storage and memristive applications. The state of 
these devices is related to the wall position and thus rapid operation will require a controllable onset of 
domain wall motion and high speed wall displacement. These processes are traditionally driven by spin 
transfer torque due to lateral injection of spin polarized current through a ferromagnetic nanostrip. 
However, this geometry is often hampered by low maximum wall velocities and/or a need for prohibitively 
high current densities. Here, using time-resolved magnetotransport measurements, we show that vertical 
injection of spin currents through a magnetic tunnel junction can drive domain walls over hundreds of 
nanometers at —500 m/s using current densities on the order of 6 MA/cm 2 . Moreover, these measurements 
provide information about the stochastic and deterministic aspects of current driven domain wall mediated 
switching. 



Via the injection of spin polarized current, spin-transfer-torques (STT) allow for the manipulation or 
switching of a ferromagnetic nanostructure's magnetization configuration without the use of an externally 
applied magnetic field 1,2 . One way for this switching to occur is through controlled domain wall motion, 
which can be driven by injecting current laterally along a nano-scale ferromagnetic strip containing one or more 
domain walls 3 6 . This ability to displace domain walls via current has driven strong research efforts into scalable, 
spintronic, domain wall-based devices such as shift registers 711 , binary memories 1214 and memristors 1517 . 

In the conventional lateral current geometry, the non-adiabatic spin torque that drives steady state domain wall 
motion 18,19 is intrinsically small 20 22 , generally limiting the domain wall speed to below 150 m/s 22 24 and requiring 
current densities above 100 MA/cm 2 , a value which can be prohibitive for applications 7 . Recent results have 
shown that the domain wall velocity can be strongly enhanced by the addition of a metallic layer with strong spin- 
orbit coupling below the ferromagnetic strip 25 27 . Nevertheless, in these latter systems the current densities 
required to move the wall at high speed are still high (—100 MA/cm 2 ) and the physical interpretation of the 
observed effects remain complicated since spin orbit coupling can also affect the domain wall structure and thus 
its dynamics 28 . 

Recently we proposed an alternative way to enhance the driving torque for domain wall motion through the use 
of vertical (or 'perpendicular') current injection through a magnetic tunnel junction consisting of a fixed magnetic 
layer which acts as a spin polarizer and a 'free' magnetic layer containing a domain wall. In this injection 
geometry, the current densities required to move the domain walls are reduced by up to two orders of magnitude 16 
and it has been predicted that wall velocities could exceed 400 m/s at current densities as low as a few MA/cm 2 29 . 

In this article, we carry out precise measurements on both stochastic and high speed deterministic domain wall 
dynamics under perpendicular current injection 16,30 35 by performing real-time measurements of domain 
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Figure 1 | Measurements of domain wall mediated resistance switching under field and vertically injected current, (a) Schematic of the MTJ stack with 
the domain wall in one of the two labeled, stable positions, (b) Typical resistance versus magnetic field curves during reversal of the free NiFe layer 
measured using a non-perturbing sense current of I sensc = 1 00 ;iA. The orange curve shows switching between the fully parallel and antiparallel states. The 
red curve corresponds to the changes in the device resistance upon cycling the domain wall back and forth between the two stable domain wall positions, 
1 and 2, near the terminations of the arc. (c) Current-induced back-and-forth displacements between positions 1 and 2 upon application of 20 ns long 
voltage pulses with / ~ 7.3 MA.crrT 2 . Insets: top view of micromagnetically simulated stable DW configurations, (d) Schematic of the setup used for the 
time resolved reflective measurement method, (e) Domain wall voltage traces (panel III) are obtained by recording the reflected voltage pulse during 
which the domain wall moves from one termination to another (panel I) and subtracting from it the reflected voltage pulse obtained when the wall is 
trapped at the second termination (panel II). (f) Single-shot voltage traces obtained under 20 ns long sub-critical current pulses (/ ~ 0.7/ c ). The parts of 
the traces occurring before the pulse and after the voltage signal has returned to zero are dashed to highlight the part of the trace corresponding to the 
waiting and propagation processes. The waiting and propagation times are labeled for the 'central' blue trace (t switc t, = t wait + t prop ~ 4 ns). 



wall-mediated magnetization switching within the free layer of arc 
shaped nano-scale MgO-based MTJ devices (Figs. la-e). These 
structures provide an ideal base architecture for data storage and 
memristive devices because the strong tunneling magnetoresis- 
tance 36,37 translates small domain wall displacements in the free layer 
into large resistance changes of the device as a whole 38 . These large 
resistance variations are also ideal for probing high speed magnetisa- 
tion dynamics 39 . 

Results 

The measured MTJ devices have the following general composition 
wherein a synthetic antiferromagnet (SAF) is used as a current polar- 
izer: SAF/MgO(1.17 nm)/NiFe(5 nm) (see methods for full device 
structure). Each device is in the shape of an arc of —115 nm width 
with horizontal terminations at each end (Fig. la). To create a trans- 
verse 40 domain wall (DW) within the NiFe layer, we apply a large in- 
plane field perpendicular to the long axis of the device. A domain wall 
is created at the apex of the arc and then, in zero field, relaxes to a 
stable position near one of the horizontal terminations which serve to 
constrain the domain wall inside the arc. Each one of these positions 
(referred to as positions 1 and 2) is well reproduced by micromag- 
netic simulations (see Figs, la and lc) and corresponds to a discrete 
resistance level due to the device's TMR (Figs, lb and lc). 

A prepared domain wall can be driven back and forth between 
positions 1 and 2, thus switching the device between the two discrete 



resistance states. This can be done either via the application of a field 
along the long axis of the device or vertical current injection in zero 
external field. Fig. lb demonstrates the field induced switching 
wherein the field is cycled between positive and negative values 
whilst monitoring the sample resistance. Notably, after displacing 
the wall to either of the two positions, it is stable at remanence in 
zero field. The application of negative (positive) current pulses can 
also be used to drive the transition 1 — > 2 (2 — > 1 ) in zero applied field. 
Fig. lc shows a series of measurements of the device's resistance 
made directly following the application of alternating positive and 
negative vertically injected current pulses with / ~ ±7.3 MA.cm 2 
and duration 20 ns. We thus demonstrate that the domain wall 
mediated switching can be carried out repeatedly using injected cur- 
rent pulses without the use of an external magnetic field. 

In the remainder of this article, we concentrate on magneto- 
transport experiments carried out to probe the dynamics of the 
domain wall during its motion between the two terminations. 
Notably, we show that the switching can be deterministic and that 
the propagation occurs at velocities on the order of 500 m/s corres- 
ponding to device switching on the order of 1 ns with negligible 
switching stochasticity. 

A schematic of the measurement circuit is given in Fig. Id. We use 
a reflection technique where we send a pulsed current through the 
MTJ and measure the reflected voltage pulse. For a domain wall 
prepared at position 1, a negative current pulse will displace it to 
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position 2 (Fig. lc). Its displacement will result in an instantaneous 
drop in the sample resistance and thus also in the magnitude of the 
reflected voltage pulse (panel I of Fig. le). A second, identical voltage 
pulse is then sent through the sample however the wall does not 
appreciably move due to the termination -induced confinement. 
The reflected voltage pulse is again recorded (panel II of Fig. le) 
and used as a reference trace which is subtracted from the first trace 
to yield a domain wall voltage signal such as that shown in panel III of 
Fig. le. We note that the device resistance is measured before and 
after both the driving and reference pulses using a non-perturbing 
100 /(A sense current in order to verify the domain wall's position at 
each step of the experiment. 

Three measured domain wall voltage traces obtained under a sin- 
gle sub-critical current density (/ ~ 0.7/ c ) are shown in Fig. If. The 
subtracted voltage signal is zero before the pulse, and then plateaus 
during a stochastic waiting or depinning time, t wajt , wherein the wall 
does not measurably move. Finally, a drop in the signal is observed, 
corresponding to the domain wall traveling to the other end of the 
device during a time t prop . The full switching time of the device is thus 



given by f s , 



t-wait + tprop- At higher current densities, as will be 



shown below, the switching time is on the contrary determined 
almost exclusively by the propagation time corresponding to an 
abrupt change in the voltage trace. 

The voltage drop corresponding to the domain wall propagation 
occurs rapidly with t prop on the order of 1 ns. To obtain a clearer 
visualization of these fast propagation dynamics, we recorded —100 
real time voltage traces, and plotted them as V versus t — t swjtch where 
tswitch was measured manually. The curves were then averaged to 
reduce the signal noise during the propagation part of the curve. 



These averaged shifted traces for three values of / are shown in 
Fig. 2a from which an upper bound for t prop can be estimated to be 
— 1.5 ns. This corresponds to an average domain wall velocity on the 
order of 400 m.s -1 given that the propagation distance is —600 nm. 
This result can be compared to those obtained from micromagnetic 
simulations (see methods) shown in Fig. 2b which show the spatially 
averaged longitudinal component of the free layer magnetization 
during current driven domain wall motion. The plotted magnetisa- 
tion component is representative of the resistive change induced 
during domain wall motion. By also shifting these simulated curves 
by t sw i tc f, it can be seen that the simulated and measured propagation 
times are in good agreement. Moreover, the simulation indicates that 
the maximum velocity during propagation can reach 600 m/s. 

Despite these fast dynamics, at sub-critical current densities the 
waiting times preceding the domain wall propagation are character- 
ized by significant variability, as seen directly in Fig. If. However, 
these waiting times become immeasurably small above the critical 
current density. To demonstrate the transition from sub-critical 
thermally activated domain wall dynamics towards deterministic 
dynamics, we measure the characteristic domain wall mediated 
switching time as a function of / as extracted from averaged voltage 
traces (—200 single shot traces per averaged trace). The averaged 
traces were obtained under 20 ns long pulses for different current 
densities with / > 0 and are shown in Fig. 2c. Here, there was no 
horizontal shift of the traces along the time axis (as done in Figs. 2a 
and 2b) and so these averaged traces probe the entire switching 
process: both stochastic depinning and propagation. At low current 
densities it is clear that the characteristic waiting time is much longer 
than the propagation time (typically 1-1.5 ns as shown earlier). 



a 








> 


5 ; 








4 - 


<D 




o> 


3 - 


re 




+-» 


2 - 


o 




> 


1 - 


"(5 




c 


0 - 


Ul 




'</5 


-1 - 


1 


-2 - 







Lprop 







" 5.7 MA. cm' 2 




6.6 MA. cm 2 




_ 8.8 MA. cm 2 













-2.0 



-1.5 



-1.0 



-0.5 



0.0 




time [shifted] (ns) 



-1.5 -1.0 -0.5 

time [shifted] (ns) 



> 
E. 

a> 
u> 
ns 

o 

> 



Cf) 

'(/> 




5 10 

time (ns) 



J2 
TO 

o 
a. 

O) 



a. 

■D 



1.0 
0.8 


o 


o 

0 

o 

o 


COOOOOOOO O O O 0 
0 ° 

o 

O o o 


0.6 


■ % 

8 


o 
o 


o o °°o 
°o 


0.4 


o 


°o° 






o 




0 -4.6MA.cm 2 


0.2 


- «? 

: $* 




0 - 5.1 MA.cm J 
- 6 MA. cm 2 


0.0 




1 1 





5 10 15 20 

pulse duration (ns) 



Figure 2 | Domain wall propagation, switching and depinning times, (a) The propagation time t prop is measured by averaging several single-shot 
traces shifted by t — t switc f,. The resulting traces are shown for different current densities and reveal t pmp — 1 ns. (b) Micromagnetic simulations: the 
average x component of the magnetization is equivalent to the resistive change due to the domain wall displacement, (c) Averaged single shot DW voltage 
traces for several current densities at a constant pulse duration of 20 ns. (d) Depinning probability as a function of pulse duration for different current 
densities in a second device. 
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Indeed the low J traces can be decomposed into two parts. First 
because of the low probability of switching events, the voltage 
remains constant. The voltage then decreases towards V = 0 with 
a smeared, sloping transition indicative of a wide distribution of t wait 
times which dominate the switching characteristics and limit switch- 
ing reliability. As the current density is increased however, the 
switching transition is faster and the waiting time reduced, resulting 
in strongly diminished plateau widths at the beginning of the aver- 
aged traces. At current densities above approximately 6 MA.cm -2 
(green and black curves in Fig. 2c) there is no measurable waiting 
time and the motion is consistent with a deterministic motion with 
tswitch ~ tp r op- Additional averaged traces are given in Supplementary 
Fig. SI. 

Conventional depinning probability measurements were carried 
out on a second device with the same structure where the probability 
of observing a depinned domain wall was measured following the 
application of a pulse with a given current density and variable dura- 
tion. Data is shown in Fig. 2d where we plot the depinning probabil- 
ity versus pulse duration. Although these measurements only probe 
twaiu we observe consistency with the the averaged voltage traces 
presented in Fig. 2c in that for a given pulse duration, the depinning 
probability increases with /. Alternatively, it is possible to normalize 
and invert the domain wall voltage traces from Fig. 2c to generate 
curves which show the probability of the domain wall having both 
depinned and traveled to the other side of the device under the action 
of a continuously applied current pulse. When tprop < < twait (ie- at 
low J), the averaged voltage traces and conventional depinning prob- 
ability measurements compare well (Fig. 3a). At higher /, the depin- 
ning probability increases faster than the switching probability 
(Supplementary Fig. S2), consistent with the fact that the switching 
measurements probe both the depinning and propagation times. 

Since it is the full switching time, tp rop + t wait , which is of import- 
ance for device operation, we use the inverted, averaged domain wall 
voltage traces to extract a mean switching time x as a function of /. To 
do this, we fit the inverted and rescaled data with a statistical cumu- 
lative distribution function which describes the probability of the 
device having switched at a given time t. The existence of a plateau 
at the start of the low / averaged traces does not permit the use of the 
exponential probability distribution function, often employed for 
magnetic switching in cases where t prop -C t wa i t 41 ' 42 . Instead, we use 
the gamma statistical distribution, a flexible 2 parameter model 
which can be employed to model waiting times 43 . Notably this dis- 
tribution can be used to fit the broad voltage traces observed at low / 
where stochasticity dominates as well as the sharply varying traces 
obtained at high / in Fig. 2c. The extracted values of x (see methods) 
are given in Fig. 3b where we plot In x versus /. A simple linear fit to 



the low / data highlights the existence of two regimes, one below and 
one above a critical current density of ~6 MA/cm 2 . This value 
appears to be quite consistent with the averaged voltage traces where 
we see that at approximately 6 MA/cm 2 the plateau in the averaged 
traces disappears. Furthermore, at this point, the drop in the voltage 
occurs over a time frame consistent with the measured and simulated 
tp rop values (Figs. 2a and 2b). Indeed, at the highest measured current 
density, the domain wall signal voltage drops to zero in a time close to 
1 ns, indicative of domain wall velocities approaching 600 m/s. A 
more rigorous analysis of the In x versus / data is discussed below 
where we extract a value for the effective field arising from the 
STT. 



Discussion 

Finite, pre-switching waiting times have been observed both below 
and above the critical current for STT induced uniform magnetic 
switching 42,44 , the latter being due to a zero torque when the magne- 
tizations of the free and polarizing layers are aligned. In our system 
however, the spins in the transverse domain wall always form an 
angle with the magnetization of the polarizing layer, resulting in a 
finite initial torque. The waiting times for domain wall propagation 
that are observed here thus simply indicate that the applied current 
density is sub-critical, meaning that thermal activation is required to 
overcome the domain wall pinning energy and initiate motion. 

The / dependence of x can be understood by considering a trans- 
ition from thermally activated to deterministic domain wall 
dynamics 42 . The characteristic time associated with thermally acti- 
vated depinning can be expressed as 

tgc exp(E B /k B T) (1) 

with the following expression for the energy barrier 41,45 

E B x(l-Hj/H de p yT=0 )K (2) 

which is defined only for Hj £ H dep T=0 . H de „ T=0 is the depinning 
field at zero temperature. Hj = H S tt + H 0e is the equivalent field 
corresponding to the net effect of the injected current, ie. the out-of- 
plane torque 29 , H STT , and the Oersted field 16 , H 0e , which both drive 
domain wall motion in our system. In the plot of x versus / in Fig. 3b 
we can identify two regimes. In the under-critical regime J < J c , the 
domain wall overcomes the energy barrier E B with the help of ther- 
mal activation resulting in a distribution of waiting times as a con- 
sequence of the stochastic depinning process. The low / data in this 
regime is well fitted using Eqs. (1) and (2) from which a value of ] c = 
5.9 MA.cm -2 can be directly extracted assuming an equivalence of / 
and H. When the current is strong enough, Hj equals H depT=0 and 
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Figure 3 | Determination of the critical current for the transition from thermally activated to deterministic domain wall motion, (a) Comparison 
between the depinning data, switching probability data and a Gamma cumulative distribution function fit (/ = —5.45 MA/cm 2 ). (b) Mean switching 
times extracted by fitting the Gamma cumulative distribution function to inverted and scaled time resolved DW signal voltage traces, (c) SEM image of 
the sample and simulated Oersted field taking into account the contribution from the current distribution within the electrodes. 
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the depinning process becomes deterministic with no appreciable or 
measurable waiting time. Even above ] c , % continues to decrease with 
increasing /, which is indicative of a faster movement out of the stable 
domain wall position (as observed in the simulations shown in 
Fig. 2b) and thus a reduced t switch ~ t prop . 

It is possible to extract the efficiency of the spin transfer out-of- 
plane torque from the thermally activated model using Hj — ^stt 
H 0e . We have calculated the amplitude and spatial distribution of 
H Qe by taking into account the actual current distribution through 
the MTJ stack and the electrodes (see Supplementary Note SI). The 
net effect of H 0e was found to aid domain wall depinning. Although 
the Oersted field can be rather strong in the terminations (up to 
60 Oe for 7=6 MA.cirT 2 , see Fig. 3c), its net effect is quite low at 
the pinned domain wall position. Simulations of domain wall depin- 
ning from position 1 under the influence of H 0e only yielded an equi- 
valent magnetic field per unit current density of 1.12 Oe.MA^'.cm 2 . 
Therefore the Oersted field results in an effective field acting on the 
domain wall of H 0e — 6.6 Oe at J c . A lower bound was also estimated 
for the zero temperature depinning field by extrapolating temper- 
ature dependent depinning measurements with H dep \ T=0 ~ 27 Oe 
(see Supplementary Note S2). As found previously in quasi-static 
measurements 16 , this value confirms the two orders of magnitude 
gain compared to domain wall motion induced by non-adiabatic 
torque in the traditional lateral injection geometry, for which the 
equivalent field at the same current density is typically 0.1 Oe 46 . 

We have experimentally demonstrated domain wall propagation 
in a MTJ at velocities on the order of 500 m.s~ 2 under perpendicular 
current injection in the absence of an applied magnetic field. Real 
time magnetotransport measurements have allowed for an analysis 
of the transition between stochastic and deterministic domain 
wall dynamics which is crucial for any type of domain wall based 
device. At current densities above the critical current density 
(~6 MA.cm~ 2 ), the device switching proceeds quickly with no mea- 
surable waiting time. This performance compares well with state-of- 
the-art uniform switching in MTJs 42 and with theory 2 " suggesting 
that a factor of about five in efficiency could be gained by using a 
polarizer which is perpendicular to the magnetization of the free 
layer. In future applications, device sizes could be further reduced 
using materials with perpendicular magnetic anisotropy to decrease 
the domain wall width by an order of magnitude. This could allow for 
device dimensions on the order of 25 X 100 nm 2 leading to sub mA 
critical currents. The vertical current injection geometry in magnetic 
tunnel junctions is therefore extremely promising for future high 
speed domain wall based spintronic devices such as memristors 
and magnetic memories. 

Methods 

Samples. MTJ devices were fabricated from a multilayer film having the following 
composition PtMn(15 nm)/CoFe(2.5 nm)/Ru(0.9 nm)/CoFeB(3 nm)/ 
MgO{1.17 nm)/NiFe(5 nm)/Ru(10 nm) which was sputtered in a Canon ANELVA 
chamber. Details of the growth and fabrication process have been presented 
elsewhere 47 . The samples were annealed at 0.5 T at 285°C during 6 hours and 
patterning carried out using electron beam lithography and ion beam etching. The 
tunneling magnetoresistance of the patterned devices is on the order of 25%, with a 
low RA product of — 12 0.41m 2 . In our convention, a positive current favors the 
anti-parallel state. 

Micromagnetic simulations. Dynamic simulations of the domain wall motion under 
the influence of perpendicularly injected current were carried out using the finite 
difference micromagnetic code, SpinPM, which was developed at Istituto PM srl. The 
simulated device had the geometry extracted from the SEM image in Fig. 3c and took 
into account the nonvertical device edges induced by the etching process. The 
following magnetic parameters were used: cell size of 3 X 3 X 5 nm 3 , a = 0.01 for the 
magnetic damping parameter, M s — 500 emu/cm 3 for the free layer magnetization 
and respectively M s = 740 emu/cm 3 and 690 emu/cm 3 for the top and bottom 
ferromagnetic layers in the synthetic antiferromagnet. The interlayer exchange 
coupling between the top layer of the synthetic antiferromagnet and the free layer was 
set to 0.0025 erg.cm -2 . The coupling between the layers within the synthetic 
antiferromagnet was set to —0.1 erg.cm -2 . The current was assumed to be uniform 
throughout the structure and a spin polarization of 0.5 was used. 



Gamma distribution. The Gamma probability distribution function is given by 



p(t)= 



'exp(-f//?) 



(3) 



where T(a) = §™ y a 1 exp( —y) dy is the conventional Gamma function, x and 
characterize the form of the distribution with the Gamma and exponential 
distribution functions being equivalent for a — 1. The inverted and normalized 
averaged voltage traces were fitted with the cumulative distribution function of the 
gamma distribution using the CDF, GammaDistribution and FindFit functions built 
into Wolfram Mathematica. This allowed a and to be determined for each value of /. 
We then calculated the mean of the distribution, aft, in order to obtain a value for t = 
a/5. A fit of one of the inverted and normalized averaged voltage traces is given in 
Fig. 3a. 
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